Magnetization and spin dynamics of a Cr-based magnetic cluster: CryNi 
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We study the magnetization and the spin dynamics of the CryNi ring-shaped magnetic cluster. 
Measurements of the magnetization at high pulsed fields and low temperature are compared to 
calculations and show that the spin Hamiltonian approach provides a good description of CryNi 
magnetic molecule. In addition, the phonon-induced relaxation dynamics of molecular observables 
has been investigated. By assuming the spin-phonon coupling to take place through the modulation 
of the local crystal fields, it is possible to evaluate the decay of fiuctuations of two generic molecular 
observables. The nuclear spin-lattice relaxation rate 1/Ti directly probes such fluctuations, and 
allows to determine the magnetoelastic coupling strength. 
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I. INTRODUCTION 



The great efforts devoted to the synthesis and investi- 
gation of nanosize magnetic molecules are motivated both 
by interests in fundamental physics and by the envisaged 
technological applications. For instance, some of these 
systems have shown phenomena such as quantum tun- 
neling of magnetization between quasi-degenerate levels, 
slow relaxation at low T, and revealed to be promising for 
high density information storage and quantum comput- 
ing [l], 0, The magnetic core of molecular magnets is 
constituted by transition metal ions sorrounded by an or- 
ganic shell which prevents intramolecular magnetic inter- 
actions. As a result, the microscopic properties of these 
nanoscale clusters can be investigated by means of bulk 
samples. Among these systems, there are homonuclear 
antiferromagnetic (AF) ring-shaped molecules formed by 
n transition metal ions in an almost planar ring. In par- 
ticular, in even membered rings the dominant AF ex- 
change interactions lead to a singlet S't = ground state 
and the energy spectrum is characterized by rotational 
bands, with the lowest-lying levels approximately follow- 
ing the so-called Lande's rule In this paper we study 
the magnetization and the phonon-induced relaxation in 
the heterometallic ring CryNi. This compound derives 
from the even membered AF ring Crg and thus provides 
a opportunity of a deeper insight in the role of topology 
in the static and dynamical quantum properties of mag- 
netic wheels [1, 0| • CryNi compound is obtained by the 
chemical substitution of a Cr'^^ ion with a Ni^+ ion in the 
structure of the Crg ring. This leads to a new molecular 
system formed by an odd number of unpaired electrons 



with dominant AF nearest neighbour interactions as in- 
ferred by susceptibility measurements 1.4]. The resulting 
S't = 1/2 ground state has been shown to be suitable to 
encode a qubit b\. 



II. MAGNETIZATION 

The magnetic molecule has been theoretically analyzed 
within a spin Hamiltonian approach, with the Hamilto- 
nian given by: 
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where s(i) is the spin operator of the the ith ion in the 
molecule (s(i)=3/2 for Cr^+ ions, and s{i)—l for the Ni^+ 
ion) . The first term of the above equation is the dominant 
nearest neighbour Heisenberg exchange interaction. The 
second and the third terms describe the uniaxial local 
crystal fields and anisotropic intracluster spin-spin inter- 
actions respectively (with the z axis assumed perpendic- 
ular to the plane of the ring). The last term represents 
the Zeeman coupling with an external field H. The pa- 
rameters of the above Hamiltonian were determined by 
inelastic neutron scattering (INS) experiments [!,[§]. In 
order to corroborate the microscopic description of the 
CryNi from INS data, a detailed study of high field mag- 
netization is very powerful. In fact, with high pulsed 
fields up to almost GOT, spin multiplets not accessible to 
a standard INS experiment can be explored. In Fig[T^ the 
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III. SPIN DYNAMICS 




Decreasing Field 

Calculation 

Increasing Field 

H ^ \ ^ \ ^ \ ^ \ 1 



(b) 




FIG. 1: (Color online) Magnetization curve of CryNi (top) 
and derivative dM/dH (bottom) at T=1.3K. Red (dark gray) 
and black lines represent the down and up experimental mag- 
netic field processes respectively. The dashed blue lines rep- 
resent the theoretical calculation with the following parame- 
ters: Jcr-Cr=16.9K, Jcr-jvi=19.6K, dcr=-0.3K, djvi=-4K, 
ffcr=1.98, gm=2.2. 



magnetization curve as a function of the magnetic field 
H is reported. A clear staircase structure with plateaus 
at ~ l/^B, ~ and odd multiples of /is reflects the 
change in the ground state due to the external field at 
the level anticrossing fields An hysteresis of the mea- 
sured magnetization curves has been observed. The ef- 
fect arises from the non-equilibrium condition due to the 
high pulsed magnetic field with a few millisecond dura- 
tion [7| and has been discussed in terms of phonon bottle- 
neck effects and magnetic Foehn effects There is 
a very good agreement between the measured and cal- 
culated magnetization curves. This is clearly visible in 
FigHb where the positions of the main peaks of the cal- 
culated and measured dM/dH matches correctly. The 
smaller peaks in the experimental dM/dH are due to 
level anticrossings between excited energy levels. The 
effects are caused by the non-equilibrium exeperimental 
conditions and are not included in equilibrium calcula- 
tions reported in Fig[T] Q. These results confirm that 
the microscopic picture derived from INS experiments 
[1, Q perfectly holds even at very high applied magnetic 
fields. 



A major obstacle to the proposed technological appli- 
cations of magnetic molecules is constituted by phonon- 
induced relaxation. In fact, molecular observables, e.g. 
the magnetization, are deeply affected by the interac- 
tion of the spins with other degrees of freedom such as 
phonons ^IQj . Here we investigate the molecular spin-spin 
correlations through an approach based on a density ma- 
trix theory . The irreversible evolution of the density 
matrix p{t) can be determined through the secular ap- 
proximation and focusing on time scales detectable by 
low-frequency techiques such as NMR. Within this the- 
oretical framework, a general expression for the quasi- 
elastic part of the Fourier transform of cross correlation 
functions is given by [l^, • 
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where N is the dimension of the Hilbert spin space of the 
molecule, Pm*' is the equilibrium population of the mth 
level and Bmm = {m\B\fn), |m) being the mth eigen- 
state of the spin Hamiltonian, while W is the so called 
rate-matrix. The mn element Wmn of W represents the 
probability per unit time of a transition between eigen- 
states |m) and \n) induced by the interaction of the spins 
with phonons. By assuming that spin-bath interaction 
takes place mostly through modulation of local crystal 
fields, the rate matrix can be calculated on the basis of 
the eigenstates of molecular spin Hamiltonian by first- 
order perturbation theory. With the choice of a spheri- 
cal magnetoelastic (ME) coupling [I^l the transition rates 
Wmn are given by: 
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X {n\Oq,,q,{sj)\m) , (3) 

with n{x) = (e'^^" - Amn = {Em - En)/h the gap 

between the eigenstates |m) and \n) of the molecule. In 
the last equation Oo,.a^(s^) = {Sqi^iSq2,i + Sg2,iSgi,i)/2 
are quadrupolar operators . Finally, 7 represents the 
spin-phonon coupling strength, which can be determined 
by comparing the theoretical results with experimental 
data. In fact, the nuclear spin-lattice relaxation rate 1 /Ti 
probes the fluctuations of molecular observables, thus 
giving information on the relaxation dynamics [10]. Ex- 
ploiting the Moriya formula [15j , the proton NMR 1 /Ti 
can be evaluated in absolute units using as inputs the 
positions of the Cr and Ni ions and of the hydrogens of 
the molecule: 

^= E <'(5,^^,,(.,) + 5^,^^,(-.,)), (4) 
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FIG. 2: (Color online) Experimental data (scatters) and cal- 
culations (lines) of reduced proton NMR l/{TixT) for differ- 
ent values of the applied field along z (parallel to the ring 
axis). 



tocorrelation function of M Ss s (w,T, H) = 

Er=i,N M>^^.T, H)A,(r, H)/[A,(T, Hf+J]. This equa- 
tion shows that the spectrum of fluctuations of M is given 
by a sum of N Lorentzians, each with characteristic fre- 
quency Xi, given by the eigenvalues of — M^. For a wide 
range of H and T in these systems only a single relaxation 
frequency Aq significantly contribute to iSs^^s^ (cj, T, H). 
As a result, if the dominant frequency Aq intersects the 
Larmor angular frequency, i.e. when Ao(To) — uj^, at 
the temperature Tq the proton NMR 1/Ti shows a sharp 
peak Being an heterometallic ring, this explana- 

tion does not hold for CryNi and Eq.(|3]) has to be used. 
Nevertheless, our calculations show that a peak in the re- 
duced l/(TixT) occurs in agreement with experimental 
data (see Figl2]). By fitting the observed peak position 
we have obtained 7 = 0.8 x 10~''THz~^. 



IV. CONCLUSIONS 



where the S , <,'{ujl) are the Fourier transforms of the 
cross correlation functions from Eq. ([2]) calculated 

at the Larmor angular frequency ojl, while the a^J 
are geometric coefficients of the hyperfine dipolar inter- 
action between magnetic ions and protons probed by 
NMR. The occurence of a peak in the proton NMR 
l/Ti has been clearly explained in homonuclear ring- 
shaped molecules with small anisotropy such as Crg 
p^ . In fact, in this case l/Ti oc Ss^.sA^l), where 
S's^^Sj (1^; H, T) is the Fourier transform of the au- 



A magnetization study of the heteronuclear antiferro- 
magnetic ring-shaped nanomagnet Cr7Ni has been per- 
formed. A clear step- wise increase of magnetization with 
increasing field is observed. The very good agreement of 
high field magnetization measurements up to almost 60T 
with calculation shows the spin Hamiltonian approach to 
be suitable even at very high fields. The relaxation dy- 
namics of the compound has been investigated by the 
proton nuclear-spin relaxation rate 1/Ti. Our calcula- 
tions are in very good quantitative agreement with ex- 
perimental data. 
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